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We give a brief sketch of lattice structure function calculations and review previous 
results for the axial coupling gA- We outline a new technique for treating fermions on the 
' lattice that preserves chiral symmetry, domain wall fermions. Finally, we give preliminary 

| results for the nucleon spectrum using this new technique. Remarkably, a large mass splitting 

between the N and N* , roughly consistent with experiment, is produced in the calculation. 
These results are encouraging for proposed calculations of nucleon structure functions. 
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§1. Introduction 

Structure functions are an important quantities for understanding hadronic 
physics from first principles and for confronting experimental measurements. In 



o 

QCD, the theory of the strong interaction, they probe the underlying make-up of 
hadrons from the fundamental quarks and gluons. A complete theoretical under- 
standing of hadronic physics requires calculations of these important functions in 
QCD (or equivalently, a set of operator matrix elements). Because QCD is a con- 
fining theory, the complete calculation of the structure functions is necessarily non- 
perturbative. A rigorous, systematic method for calculating non-perturbatively in 
QCD is the lattice method, or lattice QCD. 

The definition of nucleon structure functions follows from the factorization of 
the deep inelastic scattering (DIS) cross section of electrons and nucleons. 

£ ' °dis ~ (14) 



W ^ = j-J & A xe iq - x (P\[J»(x), J v {$j\\P) (1-2) 



where L^ v is the lepton tensor, W^ v is the hadronic tensor, and J M is the electro- 
magnetic current. We are interested in the form of W^ v 1 \ 



^-^)^-^f L]F ^ ) - (1 ' 3) 
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where F±(x, Q 2 ) and F2(x, Q 2 ) are the unpolarized nucleon structure functions. Sim- 
ilarly, for scattering off of polarized nucleons with spin S, we have 

W*(p,q,S) = i^paj^ (V 5 i(*,Q 2 ) + (V - j^j 52 (z,Q 2 )) , (1-4) 

where g\(x,Q 2 ) and g2(x,Q 2 ) are the polarized nucleon structure functions. In the 
above, Q 2 = —q 2 is the momentum transfer from the electron to the struck nucleon, 
p^ is the incoming nucleon momentum, and x = Q 2 / (2p • q). 

The structure functions F±, F2, <?i, and 52 are calculable from first principles in 
QCD. One simply needs to calculate the matrix element in Eq. 1-2. 

The nucleon state \P) is a low energy hadronic bound state of quarks and gluons; 
thus the matrix element, Eq. 1-2, cannot be calculated in perturbation theory. We 
must use a non-perturbative technique such as the lattice method to calculate it. 
However, the lattice method works only in Euclidean space-time, while DIS occurs 
near the light-cone in Minkowski space-time. Thus we cannot calculate the desired 
matrix element directly. Instead, we resort to the operator product expansion; the 
non-local operator [J M (x), J u (0)] is given in terms of an infinite set of local operators. 
The utility of the method is that most of these operators are suppressed by powers of 
1/Q. The organization of the expansion is in terms of the so-called twist of the (local) 
operator, which is defined as the dimension of the operator minus its spin, T = d — s. 
The Wilson coefficients of the expansion, C n are calculable in perturbation theory, 
and one can show from dimensional reasoning that C n ~ q~^ T ~ 2 \ so at leading order 
we need only consider T = 2. We note that the structure functions can be related to 
T = 2 distribution functions in the parton picture, i.e., the probability of finding a 
parton (quark or gluon) with momentum fraction x of the nucleon inside the nucleon. 
At tree level, the coefficient functions are simply the charges of the partons. 

In summary, the calculation of the matrix element is split into two parts, a soft 
part given by the matrix elements of the local operators, and a hard part given by 
the C n . The dividing line, or scale ji, where the split is made is arbitrary. However, 
at any finite order of perturbation theory, there is residual scheme dependence in 
C n , so one cannot take fi too small. In addition, present lattice simulations have a 
cut-off I /a around 2-3 GeV. These considerations constrain the range of fi values 
that can be used. 

Looking in more detail, we see that the matrix elements calculated on the lattice 
are actually moments (in x) of the structure functions, or more precisely, their Mellin 
transforms. Thus in principle, to construct the structure functions one must calculate 
all of the moments and perform the inverse Mellin transform. Fortunately, for x not 
very close to 1, we only need the lowest few moments since rapidly, the higher 
moments begin to probe only the high x region. 

The axial and tensor charges are given by the lowest moment of the distribution 
functions gi(x) and h±(x), the probability of finding a parton with longitudinal 
or transverse polarization inside a longitudinally or transversely polarized nucleon, 
respectively. 

Jdx gi (x)=Aq, (1-5) 
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2M N S li Aq = (P, S\q W q\P, S), (1-6) 

Jdxh 1 (x) = 5q, (1-7) 

2(^P„ - P„S„)(fg = (P, S\qia^ l5 q\P, S), (1-8) 

where Z\I7 = 1/2 J] 9=UjdjS Z\g may be interpreted as the quark contribution to the 
nucleon spin. Also, Au — Ad = g^, the axial coupling. To date, lattice calculations 
have underestimated the experimental value of the axial charge, §a = 1.26, by 
about 20-25%. The tensor charge, which is the lowest moment of the transversity 
function hi 2 \ is not yet known from experiment, but may soon be measured in 
polarized Drell-Yan experiments at Brookhaven's Relativistic Heavy Ion Collider 3 ). 
Theoretical predictions are thus desirable. So far, Kuramashi has given the only 
estimate 4 ) . 

Another detail to note is that lattice and continuum operators are obviously 
regularized in different schemes, so in any calculation that relies on both regulariza- 
tions, one must match the two, or require them to be equal, at some renormalization 
scale, ji. 

Oijj) = Z ( af i)O lat (a) (1-9) 

This matching can be done completely in perturbation theory, but it is well known 
that the bare lattice coupling constant is a poor expansion parameter and lattice 
perturbation theory is more difficult, besides. An attractive alternative is to use a 
non-perturbative renormalization (NPR) procedure 5 ). The main idea is to mimic a 
continuum momentum (MOM) scheme on the lattice. Recent calculations show this 
procedure works quite well 6 ) . 

§2. Domain wall fermions 

A promising new approach for treating fermions on the lattice is domain wall 
fermions (D WF) 7 ) . DWF maintain the full chiral symmetry of the continuum at 
finite lattice spacing. Conventional lattice fermions explicitly break chiral symmetry, 
an artifact of removing the lattice doublers, which is only restored in the continuum 
limit. By maintaining this crucial continuum symmetry, simulations with DWF are 
more continuum-like. 

The exact lattice chiral symmetry arises from the addition of an extra, infinite, 
fifth dimension. The five dimensional fermions have a mass M(s) in the shape of a 
domain wall (M(s) < for s < 0, and M(s) > for s > where s is the coordinate 
in the fifth dimension). Four dimensional chiral zero modes naturally arise on the 
defect where M(s) = 0. 

The chiral zero modes are four dimensional in the sense that they propagate 
along the domain wall in ordinary 4d space-time, but not into the fifth dimension. 
In general they have an exponentially decreasing wave function in the fifth dimension. 

For a periodic extra dimension, an anti-domain wall also appears. The chiral 
zero modes on the anti-domain wall have the opposite handedness as the ones on the 
domain wall. By coupling both modes to the same 4d gauge field, we can construct 
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a vector gauge theory, for example QCD. When the gauge fields are explicitly four 
dimensional, one can think of the extra dimension as an internal flavor space. This 
idea led to the overlap formulation of DWF 8 ) . 

As long as the left and right handed zero modes do not overlap in the middle of 
the fifth dimension, the chiral symmetry is exact. On a finite lattice, however, the 
tails overlap, and the left and right handed modes mix; an intrinsic quark mass is 
generated. This quark mass can be made arbitrarily small by increasing the size of 
the 5th dimension, N s (the number of sites in the fifth dimension). 

In practice, we use the boundary fermion variant of DWF 9 '. Here, half of 
the fifth dimension is discarded, so the domain walls become the boundaries of the 
fifth dimension. An explicit quark mass is introduced by coupling the boundaries 
with a strength —m, i.e., the chiral projections of the five dimensional fields on the 
boundaries couple to each other with strength m, just like a usual fermion mass 
term. As long the intrinsic mass term is much smaller than m, it can be neglected. 

It turns out that the exponential rate of damping of the zero mode wave function 
is controlled by the 5 dimensional quark mass, or domain wall height M. Thus, there 
are three parameters controlling explicit chiral symmetry breaking, m, M, and N s . 
The formal chiral limit is N s — > oo and m — > 0. 

Because they are protected by chiral symmetry, DWF have only 0(a 2 ) errors 
io),n)_ xhis is true at any coupling and in the N s — > oo limit. At finite N s this is 
expected to hold up to exponentially small corrections. 

In the last section we will see an example of the power of DWF in QCD simu- 
lations, the calculation of the negative-parity- nucleon/nucleon mass splitting. This 
calculation also shows that DWF fermions may be a powerful tool for structure 
function calculations. 

§3. Review of lattice results for g^. 

In this section we briefly review the current results for gA from lattice calcula- 
tions. This is perhaps the simplest nucleon matrix element one can imagine, and it 
is troubling that lattice predictions are too low (see table I). Naively, one expects 
quenching effects to be small since only the valence quarks contribute in the isospin 
limit. In fact, recent dynamical results also significantly under predict gA 12 ^- 

The problem may be related to the renormalization of the lattice axial current 
Za, which in general is not equal to one because of explicit chiral symmetry breaking. 
The values of g A in table I were obtained using Wilson fermions and perturbative 
calculations of Za ■ These simulations are far from the continuum limit where lattice 
perturbation theory is trustworthy. There may also be significant discretization 
errors since each of the studies were done at a single, rather large, lattice spacing. 

For DWF, the conserved axial current 16 ) receives no renormalization. This is 
not true for the local current. However, recent calculations 6 ) show that Za/Zv = 1, 
a condition required by the axial Ward-Takahashi identity. Thus, a reliable estimate 
of g A using DWF is given by the ratio of the (unrenormalized) matrix elements 
(P\A,i\P) I \P\Vp\P) since {P\V^\P) = 1 in the continuum. 
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Table I. Summary of renormalized results for g A . Results are for Wilson fermions. In cases where 
point-split(PS) currents were used, the renormalization was calculated using local currents. 
While the second entry is consistent with the measured value of g A = 1.26, the configuration 
ensemble is quite small. Subscript "con" refers to the connected contribution to the matrix 
element. 
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§4. Nucleon spectrum using domain wall fermions 

In Fig. 1 we show the low energy nucleon spectrum as a function of the quark 
mass, to 17 ). The data are all for 6/g 2 = 6.0 (quenched), M = 1.8, and lattice size 
16 3 x 32 x 16. We use 205 gauge configurations for the lightest two quark masses, 
to = 0.02 and 0.03 and 24 configurations for the heavier ones, to = 0.04 — ► 0.125. For 
an explanation of the baryon interpolating operators Bf 2 , see Ref. 17 \ We omit the 
point at m=0.02 for the operator B 2 since a good plateau in the effective mass plot 
is absent. The experimental values for iV(939), iV'(1440), and iV*(1535) (to ~ 0) 
(bursts) have been converted to lattice units (a -1 wl.9 GeV from am p = 0.400(8) 
in the chiral limit) 18 ) . 

The most remarkable feature in Fig. 1 is the splitting between the nucleon iV 
(from the Bf operator) and its parity partner, N* (from both B^[ and B 2 operators). 
In the chiral limit, the mass difference is roughly consistent with experiment (within 
15%). To our knowledge, this is the first such demonstration of this important feature 
of the low-lying nucleon spectrum. Because of the good chiral symmetry of DWF, 
this result suggests that the mass splitting is due to spontaneous chiral symmetry 
breaking. 

A remaining puzzle is that we cannot extract the nucleon mass from the uncon- 
ventional nucleon operator B 2 which does not have a non-relativistic limit. Instead, 
B 2 gives a signal for the positive parity excited state of the nucleon. To our knowl- 
edge, no other groups have succeeded in extracting a clear signal for either the 
nucleon or its excited state using B 2 . We suspect that when using DWF the op- 
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Fig. 1. The low- lying baryon spectrum in lattice units. Note the large mass splitting of N and N* 
which is within 15% (in the chiral limit) of the experimental value (bursts). The mass of the 
positive parity excited state is too high, however. 



erator B^ has negligible overlap with the nucleon itself, contrary to Wilson quarks 
where the operators B^ and B^ mix due to explicit chiral symmetry breaking 19 ). 
We confirm such a possibility by comparing the value of the TV' mass extracted above 
to the value from a two state fit to the Bf correlation function. Finally, we see no 
signal for the nucleon state in the mixed correlation function (Bf E>2 + B% ^t) which 
leads us to conclude (0\B£\N) ~ 0. 

In Fig. 2 we show the mass ratio of the baryon parity partners (win* /ttin) versus 
the pseudoscalar and vector meson mass ratio (m n /m p ). Experimental points are 
marked with bursts and correspond to non-strange (left) and strange (right) sectors. 
In the strange sector we use U and 17(1750) as the baryon parity partners and K 
and K* for the mesons. The baryon mass ratio grows with decreasing meson mass 
ratio. A naive extrapolation is quite consistent with the experimental values. For 
comparison, we show the recent result of Lee and Leinweber 2 °) which was obtained 
using an improved Wilson quark action, Dx34, on coarse lattices. 
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Fig. 2. Mass ratios of the low lying baryon spectrum. Experimental measurements are denoted by 
bursts. The improved Wilson data (diamonds and triangles) are from Ref. 20 - ) . 



§5. Summary 

The DWF results presented here are very encouraging. The negligible mixing 
between the different operators -Bi 5 2> and the good agreement with the observed 
mass splitting of the N and N* states indicate that DWF will work well for the more 
challenging calculation of nucleon structure functions. The above also confirms our 
intuition that maintaining chiral symmetry produces more continuum-like lattice 
simulations. Systematic effects due to finite volume and lattice spacing will be 
addressed in future calculations. A calculation of nucleon matrix elements using 
DWF is to be started soon. 
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